Study on the stability analysis of rainfall slope based on G-A model considering moisture content

In this paper, the moisture content on the wetting front is obtained by using the Van Genuchten (VG) model of unsaturated soil, and then the distribution of moisture content in the upper part of the wetting front is simplified as a trapezoid. The Green-Ampt (G-A) infiltration model of infinite slope with unsaturated characteristics is derived. The analytical expression of safety coefficient (FOS) of infinite slope with rainfall is solved by combining the limit equilibrium method with the unsaturated soil shear strength theory phase. The results show that: 1) compared with the traditional G-A model and the combined rectangular and 1/4 oval model, the upper part of the wetting front is simplified to a trapezoidal model, which has great advantages in infiltration rate and cumulative infiltration, especially when the slope is large or the rain intensity is heavy; 2) since the distribution of soil moisture content above the wetting front is considered, the matrix suction at the wetting front is not neglected, and the safety coefficient calculated by the method proposed in this paper is closer to the actual situation than the traditional G-A model.


Stability analysis of unsaturated soil slope
Assumptions. Due to the complexity of realistic slope conditions, to facilitate the solution, the following assumptions are made before the theoretical derivation: (1) The rainfall is uniform (i.e., the rainfall intensity is constant); (2) The soil is homogeneous; (3) Unsaturated soil with initial moisture content below the wetting front; (4) When rainfall infiltrates, there is a prominent wetting front in the slope; (5) The slope is simplified as an infinite slope. Green-Ampt model wetting front depth calculation. When the rainfall intensity is less than the infiltration rate, according to the Green-Ampt infiltration model, the infiltration rate i in the slope determined by the rainfall intensity is expressed as follows: where i represents infiltration rate, q r represents rainfall intensity, and α represents slope angle.
The relationship between rainfall duration t and infiltration depth z f can be expressed as follows: where θ s and θ i represent the saturated moisture content and the initial moisture content.
In the traditional G-A model, the infiltration rate i is expressed as follows: where K s is the saturated hydraulic conductivity, h 0 is the depth of ponding, h f represents the suction head of the wetting front, and z f represents the wetting front depth. For an infinite homogeneous slope, in the stable infiltration process, the following relations exist: Combining Eqs. (3) and (4), the relationship between wetting front and infiltration time is as follows: Therefore, when the rainfall intensity is greater than the infiltration rate, in the traditional G-A model, the relationship between the depth of the wetting front z f and the rainfall duration t is expressed as: Improved Green-A model considering moisture content. It is difficult to reach a fully saturated state during slope infiltration. The moisture content in the soil mainly depends on the soil matrix suction, in order to obtain the distribution of moisture content in the wetting area during the rainfall infiltration process, it is necessary to analyze the distribution of the matrix suction during the rainfall infiltration process. The schematic diagram of the calculation of rainfall-induced slope infiltration is shown in Fig. 1. The downward direction of the vertical slope is selected as the positive direction. To obtain the variation law of the matric suction in the wetting area, point A in the slope, point B in the wetting area and point C in the wetting front are selected as the calculation points. The driving force of soil moisture movement is mainly determined by the potential energy difference between the soil potential energy at a certain position and the reference object 19 . According to the potential energy equation, the total potential energy is expressed as follows: where K(θ) is an unsaturated permeability coefficient. The infiltration rate between A and B and between B and C are respectively: From the Eqs. (8) (9) (10) (12) and (13), it can be known that the matrix suction at point B is expressed as follows: Matrix suction at any depth in the soil infiltration profile can be calculated from Eq. (14). The Van Genuchten model is a model that describes the relationship between soil water-holding capacity and matrix suction, and this model has been recognized and widely used by scholars 20 . The specific expression is as follows: where θ(z) represents the moisture content of the soil when the depth is z, θ s represents the saturated moisture content of soil volume, θ r is the residual moisture content, h is the pressure water head, α n and m are fitting parameters.
Combining Eqs. (14) and (15), the moisture content θ(z) at any depth in the wetting area can be obtained as: The cumulative infiltration amount F of rainfall can be calculated as follows: www.nature.com/scientificreports/ Combining Eqs. (16) and (17), the cumulative infiltration amount can be expressed as follows: In order to facilitate the calculation, the above equation should be simplified. According to the Eq. (18) and referring to the results of trapezoidal moisture content distribution and rectangular and 1/4 elliptical moisture content distribution proposed by Zhang, et al. 18 , the soil moisture content above of wetting front is simplified to trapezoid, as shown in Fig. 2. The moisture content distribution of the slope can be calculated as follows: Put z = z f into the Eq. (18), the moisture content (θ zf ) of the wetting front can be obtained as follows: Simplified moisture content in the upper of the wetting front according to trapezoid, the soil moisture content θ(z) at any depth can be simplified as: Combining Eqs. (17) and (20), the cumulative infiltration is expressed as follows: The rate of infiltration i can be derived from Eq. (21): In the beginning of the rainfall, the rainfall intensity is less than the infiltration rate, the rainfall intensity controls the infiltration rate of the soil, and the infiltration rate q w is: www.nature.com/scientificreports/ As the rainfall progresses, the slope's surface tends to be ponding, and the infiltration capacity of the slope gradually decreases. According to the G-A model, the infiltration rate of the slope can be expressed as: Due to the continuous change of infiltration rate, there is a critical time t p , and the corresponding infiltration depth is z p . At this time, the component of rainfall intensity in the direction of the slope is equal to the infiltration capacity of the soil, based on the infiltration rate, the following equation is established: The infiltration depth is expressed as follows: Combined with Eqs. (24) and (26), the critical cumulative infiltration I p can be obtained: In the joint Eqs. (27) and (23), the critical time point t p can be expressed as follows: When the rainfall intensity is less than the infiltration rate of the soil itself, that is, when t < t p : When the rainfall intensity is greater than the infiltration rate of the soil, that is, t > t p , the infiltration rate of the soil is determined by the soil itself: To sum up, the relationship between the development rate of the wetting front depth and the rainfall time in the whole process of rainfall infiltration is as follows: Integrate Eq. (33) and substitute the initial boundary conditions t = 0, z f = 0, the infiltration critical conditions t = t p , z f = z p , the relationship between time and infiltration depth in the process of infiltration can be obtained: The dynamic change with time of wetting front depth during rainfall infiltration is expressed as Eq. (34). Compared with the improved G-A model of Eq. (6). The model in this paper considers the change of moisture content in the wetting area and the depth during the rainfall infiltration process. www.nature.com/scientificreports/ Calculation of stability coefficient of infinite slope. The slope infiltration is simplified as shown in Fig. 3 and a representative soil block is taken for force analysis. In rainfall infiltration, the soil weight changes with the change of moisture content. The weight of soil at any depth is expressed as follows: where γ(z) is the soil weight at the depth z of the wetting area; γ d is the dry weight of the soil; γ i is the initial weight of the soil. By combining Eqs. (35) and (20), the soil weight at any depth above the wetting front can be expressed as follows: The weight G of the representative soil strip is expressed as follows: Taking the area of L × 1 for analysis, the sliding surface shear stress τ s and normal stress σ can be obtained from the mechanical equilibrium equation as: Based on the unified effective stress principle of unsaturated soil 21 , the effective stress is expressed as follows: where σ ′ is the effective stress, u w is the pore water pressure. According to the Mohr-Coulomb failure criterion, the shear strength of the sliding surface can be expressed as follows: www.nature.com/scientificreports/ τ f is the shear strength, c' is the effective cohesion, and φ' is the effective internal friction angle. In the limit equilibrium state, the stability coefficient (FOS) of the slope is defined as the ratio of the sliding surface shear strength τ f to the shear stress τ s , expressed as follows: FOS can be deformed as follows:

Calculation example verification
The numerical simulation results of Li, et al. 22 are used to compare. The infiltration depth and FOS for the five rainfall durations of 6 h, 7 h, 8 h, 9 h, and 10 h are given by numerical simulation, where the rainfall intensity q r = 4.333 × 10 -4 m/min, the slope angle α = 30°, and the specific soil parameters are shown in Table 1.
In the figure, TGAM represents the traditional Green-Ampt model and IGAM represents the improved G-A model in this paper. It can be seen from Fig. 4 that in the early stage of rainfall infiltration, the boundary conditions of infiltration are controlled by the rainfall intensity, and the wetting front depth increases linearly with the rainfall time. The infiltration boundary of the soil is controlled by the soil's infiltration capacity in the middle and later stages of rainfall. The infiltration rate gradually decreases and finally tends to be stable, and the depth of the wetting front shows a certain degree of nonlinearity with the rainfall time. With the infiltration time increasing, the infiltration model that does not consider the water distribution above the wetting front results in slower infiltration than the actual infiltration situation. The method proposed in this paper is closer to the actual infiltration situation. This is because it is assumed that the upper part of the wetting front is saturated, the moisture content of the soil is larger, and the wetting front reaches the same depth, which requires more infiltration time, and under the same rainfall conditions, longer rainfall infiltration times are required. Compared with the G-A model that does not consider the unsaturated region above the wetting front, the method proposed in this paper is more in line with the actual infiltration situation.
It can be seen from Figs. 5 and 6 that the FOS calculated by the two models show consistent regularity with time or infiltration depth. In the early stage of rainfall, the FOS calculated by the two models both decreased rapidly with the rainfall time. The rate of decrease of the stability coefficient calculated by the two models slows down with the continuous rainfall. Compared with the moisture content distribution without considering the wetting area, the FOS calculated by the method proposed in this paper is larger in the same rainfall infiltration time. This is because when the moisture content above the wetting front is considered to be the saturated moisture content. When calculating the FOS at the wetting front, the effect of matrix suction is often ignored. Therefore, the calculated FOS is often smaller than the actual value.
The moisture content distribution of soil section. The moisture content distribution in the soil section under different rainfall times are shown in Fig. 7. When the infiltration depths were respectively 0.5 m, 1.0 m, 1.5 m, 2.0 m, and 2.5 m, the required time was respectively 2.14 h, 4.52 h, 7.00 h, 9.51 h, and 12.05 h; the time needed for the wetting front to increase the depth of 0.5 m is 2.14 h, 2.38 h, 2.42 h, 2.51 h, and 2.54 h, respectively. It can be seen from the above infiltration law that as the depth of infiltration increases, the time required for infiltration increases; that is, the infiltration rate tends to decline. www.nature.com/scientificreports/   www.nature.com/scientificreports/ Influence of different rainfall intensities. The parameters of the soil are given in Table 1, the rainfall intensity q r = 1.139 × 10 -4 m/min is selected as the low-intensity rainfall condition, and the rainfall intensity q r = 4.333 × 10 -4 m/min is selected as the heavy rainfall condition. The variation of the wetting front with rainfall time and the variation of FOS with time under two rainfall conditions were calculated respectively. The results are shown in Figs. 8 and 9: It can be seen from Fig. 8 that there is a great correlation between FOS and rainfall intensity. In the same rainfall duration, the decline of FOS under heavy rainfall is significantly faster than that under light rainfall. When the moisture content in the upper part of the wetting front is considered as a trapezoidal distribution, the decline of FOS is significantly lower than that of the traditional G-A model. It is worth noting that when the rainfall duration is more than 21 h. The FOS calculated by the traditional G-A model under light rainfall conditions is smaller than that of the assumption that the upper part of the wetting front is trapezoidal moisture content under heavy rainfall conditions. It can be seen from Fig. 9 that the infiltration considering the water distribution in the upper part of the wetting front is faster than that of the traditional G-A model, and this phenomenon is more obvious when the rainfall intensity is larger. Therefore, in the same rainfall duration, the moisture content in the upper part of the wetting front is equivalent to a trapezoidal distribution, although the infiltration is deeper, the calculated FOS is greater.
Influence of different slopes. In order to study the variation law of slope angle on infiltration and slope stability. The slope angles of 30° and 50° were selected respectively, and the variation of infiltration depth and FOS with rainfall time was calculated. The results are shown in Figs. 10 and 11.
It can be seen from Fig. 10 that with the increase of the slope angle, the wetting front depth decreases in the same rainfall duration. This is because the angle between the rain and the slope decreases as the slope angle increases, the infiltration in the direction perpendicular to the slope surface is reduced making infiltration more difficult. It can be seen from Fig. 11 that the larger the slope angle, the worse the slope stability during the same rainfall time. It is worth noting that the upper part of the equivalent wetting front is equivalent to trapezoid and the traditional G-A model method the FOS calculated is more different, when the slope changes from 30° to 50°. It can be seen from Eq. (42) that when the soil parameters are determined, the FOS of the slope only depends on the depth and slope angle of the wetting front, and the increase of the slope angle will reduce the depth of the wetting front and the slope stability at the same time. It shows that the slope angle has a significant effect on the strength of the slope during the whole process of rainfall infiltration.
Variable-rate rainfall intensity. It is obvious that the rainfall intensity cannot be constant. The analysis model of variable rain intensity was proposed by Chu through improved G-A model 11 . Assuming that R i-1 R i is  (52) , t ≥ t x Figure 11. Variation of FOS with rainfall time at different angles. www.nature.com/scientificreports/ Case 3, there is ponding at the beginning and ponding at the end, and the whole infiltration process is controlled by the slope itself and has nothing to do with the intensity of rainfall. The judgment equation is as follows: The relationship between infiltration depth z f and time t at this stage is as follows.
Case 4, when there is ponding in t n-1 , there is no ponding in t n , and the rainfall intensity of this process is less than the infiltration rate. The judgment equation is as follows: When the rainfall time is short, net infiltration will no longer occur in this process. The infiltration in this process is ignored.

Discussion
In the actual infiltration situation, during the rainfall infiltration process, the soil will present three parts (i.e., the saturated area transition area and natural moisture content area). The cumulative infiltration and infiltration rate of 1/4 oval and rectangular combined model and trapezoidal model were compared by Zhang, et al. 18 experimentally with the real infiltration. The results are shown in Fig. 12. The result shows that at the beginning of the rainfall, the equivalent trapezoid results are better than the traditional G-A model but slightly worse than the 1/4 ellipse and rectangle combination model. As the duration of the rainfall increases, the trapezoid equivalent is much better than that of combination model of rectangle and 1/4 ellipse and traditional G-A infiltration model. This is because as the duration of the rainfall increases, the soil moisture content at the equivalent 1/4 ellipse is much larger than the actual soil moisture content. However, in the trapezoidal equivalent moisture content model, the bottom of the wetting front is higher than the actual moisture content and the top of the wetting front is lower than the actual moisture content, the moisture content of top of the wetting front and the moisture content of bottom of the wetting front cancel each other out to a certain extent, making the trapezoidal equivalent total moisture content more in line with the actual situation. For the estimation of agricultural irrigation depths with short infiltration times, the rectangle and 1/4 trapezoid combination model is preferable.
, t x ≤ t ≤ t n (55) www.nature.com/scientificreports/ However, the slope calculation FOS is generally long in the selection of rainfall time. Therefore, the trapezoidal equivalent method is more suitable for slope FOS calculation. It is worth noting that, although the trapezoidal simplified model of moisture content in the upper part of the wetting front proposed in this paper is more accurate in accumulative infiltration than the 1/4 ellipse and rectangle combination model and the traditional G-A model. The cumulative infiltration is shown in Fig. 13. But compared with the real infiltration situation, the equivalent method proposed in this paper still has the disadvantage of higher cumulative infiltration. Due to the real infiltration is very complicated. The compression of the gas in the soil will produce an air entrapment effect 1 , the compaction of the soil will become denser with the depth, and the expansive soil particles in the soil will expand and compress the soil pores under the action of rainwater 23 , the wetting front is not simple, as there is the finger instability phenomenon 24 . All these will lead to the deviation between the cumulative infiltration in the actual situation and the actual situation. The effects of these aspects on infiltration have been considered. However, how to combine these effectively to get a more accurate infiltration model is the next step should be considered.

Conclusion
In this paper, the moisture content at the wetting front is calculated, based on Darcy's law and the V-G model of unsaturated soil, then the moisture content distribution in the upper part of the wetting front is simplified to a trapezoid. Combined with the slope G-A infiltration model, a rainfall infiltration model considering the change of moisture content in the wetting area of the slope with depth is deduced. Combined with the G-A equation of infinite rainfall-induced slope infiltration and the shear strength theory of unsaturated soil, the calculation formula of the FOS of rainfall-induced landslide considering the distribution of soil water above the wetting front is proposed. And applying it to the stability analysis of infinite slopes under rainfall conditions. The following conclusions are obtained: (1) Through the analysis of the analytical expression form, the analytical expression of the FOS of rainfallinduced landslide considering the water distribution in the upper part of the wetting front is mainly composed of three parts. The factors affecting the stability of the slope under the action of rainfall can be summarized as slope angle, soil mechanics parameters and effects of matrix suction in unsaturated soils. (2) Compared with the saturated moisture content and combined model of 1/4 ellipse and trapezoidal distribution in the upper part of the wetting front. When the soil moisture content in the upper part of the wetting front is considered to be trapezoidal, the rainfall infiltration of the wetting front on the slope is faster, but the cumulative infiltration is smaller. www.nature.com/scientificreports/ (3) When the moisture distribution above the wetting front is considered to be trapezoidal, the calculated FOS of the slope increases significantly, especially when the slope angle is large and the heavy rainfall